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£C) Abstract 

,__, In the paper we consider the linear underdetermined system of a 

^J special type. Systems of this type appear in non-homogeneous network 

^J flow programming problems in the form of systems of constraints and 

fH can be characterized as systems with a large sparse submatrix repre- 

senting the embedded network structure. We develop a direct method 
for finding solutions of the system. The algorithm is based on the 
i—i theoretic-graph specificities for the structure of the support and prop- 

erties of the basis of a solution space of a homogeneous system. One 
J> of the key steps is decomposition of the system. A simple example is 

ON regarded at the end of the paper. 
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1 Introduction 

The work on this paper was motivated, mainly, by the analysis of prob- 
lems of non- homogeneous network flow optimization on large data files [l]-[3], 
[5]- [7]. Our main goal was to develop an effective (direct) method for solving 
large sparse systems of linear equations with embedded network structure, 
which appear naturally, e.g. as systems of constraints, in a broad class of 
non-homogeneous network flow programming problems. 



The 'network nature' of the regarded system allows keeping data in the 
matrix-free form in the computer memory. The formulae, derived within 
the paper, are written in the component (network) form to provide clear 
approaches towards developing computational algorithms using efficient data 
structures for graph representation [1]. 

The general idea of the method is based on the following key steps: 

• Distinguishing between the network part of the system and the addi- 
tional part. The network part of the system represents a network structure 
and corresponds to the network part of the system of main constraints of a 
non-homogeneous network flow programming problem [1], and is given, tra- 
ditionally, by balance equations, written for the nodes of a network. The 
additional part of the system corresponds to the additional part of the sys- 
tem of main constraints and can have a general form. We start the solution 
by considering the network part of the system only. 

• Introduction of the support of the network for a system. The term 
'support of the network' (also referred to as network support, or support) 
is borrowed from optimization theory [2], [3] and is used here for further 
compatibility with applications in problems of non-homogeneous network 
flow programming. The actual meaning in this paper is - a set of indices of 
variables (or, in the network terms, - a set of arcs) corresponding to columns, 
which form a basis minor of the matrix of a system. We study the support 
for the network part of the system, finding the correspondence between the 
columns of a basis minor and a family of spanning trees. 

• Construction of a general solution for the network part of the system. 
We compute a basis of a solution space of the corresponding homogeneous 
system and interpret the basis vectors as characteristic vectors, entailed by 
non-support arcs. A simple approach for finding a partial solution of the 
(non-homogeneous) system is provided. 

• Decomposition of the system. We perform column decomposition of 
the system by separating the variables according to the sets - Ut, Uc an d 
Un, which consist of the arcs of the support for the network part of the 
system, cyclic arcs and non-support /non-cyclic arcs respectively; and, finally, 
sequentially express the unknowns corresponding to the sets Uc and Ut in 
terms of the independent variables corresponding to the set Un. 

1.1 General form of the system 

Let S = (I, U) be a finite oriented connected network without multiple 
arcs and loops, where J is a set of nodes and U is a set of arcs defined 
on J x I(\I\ < oo, \U\ < oo) . Let K (\K\ < oo) be a set of different 
products (types of flow) transported through the network S. For definiteness, 



we assume the set K = {1, . . . , \K\}. Let us denote a connected network 
corresponding to a certain type of flow k E K with S k = (I k ,U k ),I k C 
I,U k = {(i,j) k : (i,j) E U k },U k C U - a set of arcs of the network S 
carrying the flow of type k. Also, we define sets K{i) = {k E K : i E I k } 
and K(i,j) = {k E K : (i,j) k E U k } of types of flow transported through a 
node i El and an arc (i,j) E U respectively. 

Let us introduce a subset U of the set U, and let K (i,j) C K(i,j), (i,j) E 
U be an arbitrary subset of K(i,j) such that \K (i, j)\ > 1. 

Finally, the initial network S = (I, U) may be considered as a union of 
\K\ networks S k , combined under additional constraints of a general kind. 

Consider the following linear underdetermined system 



Y 4 - Y 4 = a ^ *£i\kEK, 


(1) 


ie/+(c/ fc ) jei~(U k ) 




Y Y A 5H- = a P' p = ^9' 


(2) 


(i,j)EU kEK(i,j) 




V x£ = zy, (i,j)eU , 


(3) 



keK (i,j) 

where /+(£/ fc ) = {j e /fc : (hJ) k e ^*} > A r (^ fc ) = G /fc : &*)* G ^ 
af , AjJ, Q!p, % G R - parameters of the system; x = (x k j, (i,j) k E U k 1 k E K)- 
vector of unknowns. 

The matrix of system ([I]) - pi) has the following block structure: 



.4 



M 

Q 

T 



(4) 



Here M is a sparse submatrix with a block-diagonal structure of size 
y \I | x 2, 1^ I such that each block represents a \I k \ x |?7 fe | incidence 

keK keK 

matrix of the network S k = (I k , U k ), k E K, namely, 

M = Mi M 2 • • • M\ K \ , where M k ,k = l,..., \K\ are blocks of ma- 
trix M; Q is a q x N^ |£^ fc | submatrix (dense, in the general case) with 

keK 

elements AjJ, (i,j) E U,k E K(i,j) ,p — l,g ; T is a [ C/o | x /, |£^ fc | subma- 

keK 

trix consisting of zeros and ones, where all the nonzero elements appear in 
columns corresponding to arcs (i,j) k ,(i,j) E U ,k E K (i,j). We assume 

that Y, l J 1 + 1 + \ U o\ < Y l^l" 

keK keK 



2 Network part of the system 

We start the solution of system (1) - (3) by considering the network part 
of the system. 

Definition 1 We call system (1) the network part of the system (l)-(3). 
Systems (2) and (3) are called the additional part of the system (l)-(3). 

Before we proceed, let us recall the following necessary and sufficient 
condition of consistency for system (1) implied by Kronecker-Capelli theorem: 

J2 a\ = 0, k G K. 

Theorem 1. (Rank theorem). The rank of the matrix of system (HI) 
for the network S = (I, U) equals \_. \I I — 1-^1- 

k&K 

Proof. Since matrix M of the system (1) has the form 

M = Mi M 2 • • • M\ K \ , where M k is a diagonal block of matrix M, k = 

\K\ 

1,...,\K\ and rank M k = \I k \-l [1] then rank M = ^rank M k = ^(|/ fc |- 

fe=l k&K 

l) = J2\I k \-\K\. □ 

k&K 

Remark 1 We assume, without loss of generality, that the rank of the 
system Q - @ 

is 2_, \I k \ — l-^l + Q + l^o| i where q + \Uq\ is a number of equations in the 

additional part (2) - (3). 

Since the matrix of system (fTl) has the block-diagonal structure, we split 
the solution of the system into \K\ solutions of (independent) systems , each 
of which corresponds to a separate block, i.e. to a fixed k € K , and has the 
following form: 

J2 x %~ H x % = a ^ ieJk ■ ( 5 ) 

jei+(u") j£ir(u k ) 

2.1 Support Criterion 

Let's define a support of the network S = (I, U) for system (fTl). 
Definition 2 The support of the network S = (I, U) for system pi) is a 
set of arcs Ut = {U? C U k , k G K} such that the system 

E 4- E 4 = °' ^i k ,keK (6) 



has only a trivial solution for U k = U^, but has a non-trivial solution for 
U k = U k ,keK\ k ; U k ° = U k ° \J(i, j) fc °, {i,j) k ° U$>, k G K. 

Theorem 2. (Network Support Criterion). The set 

Ut = {U T , k G K} is a support of the network S = (I, U) for system Q 
iff for each k G K the set of arcs U T is a spanning tree for the network 
S k = (I k ,U k ). 

Proof. Follows directly from the proof [3] for the case when \K\ = 1 and 
the block-diagonal structure of the matrix of the system (Jl]). □ 

2.2 Basis of a solution space of a homogeneous system. 
Characteristic vectors 

Before introducing the definition of a characteristic vector, let's analyze 
the structure of a network obtained by appending an arbitrary arc (r, p) k G 
U k \U T , where k G K is fixed, to the support Ut- 

For a fixed k G K we consider a network S k = (J fc , U T |J(r, p) k ), (r, p) fc G 
U k \U T , where the set U T is a spanning tree of the network S k . Appending 
an arc (r, p) k G U k \U T to the tree entails a unique cycle. We denote this 
cycle with L k p . The set Z k = {L k p , (r, p) k <E U k \ U T } is the fundamental set 
of cycles with respect to the spanning of the network S k [1]. 

Let's consider a cycle L k p , entailed by an arc (r, p) k G U k \ U T . We de- 
fine the detour direction within the cycle L k p corresponding to the arc (r, p) k . 

Definition 3 We call an arc (i,j) k G L k p , where k G K is fixed, a 
forward arc of the cycle L k p , if the direction of the arc (i,j) k is the same as 
the direction of the arc (r, p) k within the cycle L k p . Similarly, we call an arc 
(i,j) k G L k p , where k G K is fixed, a backward arc of the cycle L k p , if the 
direction of the arc (i,j) k is opposite to the direction of the arc (r,p) k within 
the cycle L k p . 

We denote the sign of an arc (i, j) k within a cycle L k p by sign(i,j) Lr p , 

r i, (ij) k EL k r + 

szgn(z }3 ) L " P = \ -l,(i,j) k EL k T -, (7) 

( Q,(i,j) k ?L k Tp 

where L k + and L k p are the sets of forward and backward arcs of the cycle 
L k with a direction corresponding to the arc (r,p) k . 

Let us give a constructive definition of a characteristic vector, entailed by 
an arc. 



Definition 4 Characteristic vector, entailed by an arc {r,p) k G U k \ U^ 
with respect to the spanning tree Uj, , is a vector 5 k (r, p) = (<&(r, p), {i,j) k G 
U k ), where k G K is fixed, constructed according to the following rules: 

• Add an arc (r, p) k G U k \ U k , to the set Uj., k G K , which is a spanning 
tree for the network S k = (I k , U k ); and thus create a unique cycle L k p . 

• Let the arc (r, p) k set the detour direction within the cycle L k p and 

S k rp (r,p) = l. 

• For cycle's forward arcs, let SfAr, p) = 1. 

• For cycle's backward arcs, let 5 k j(r,p) — — 1. 

• Let5 k J (r,p) = 0,(t,j) k eU k \L k p . 

For briefness, further in this paper, we will call a characteristic vector 
5 k (r,p), entailed by an arc (r,p) k , with respect to the spanning tree U*, a 
characteristic vector S k (r,p), entailed by an arc (r,p) k , or, simply, a charac- 
teristic vector 5 k (r,p). 

The next two lemmas state the essential properties of characteristic vec- 
tors. 

Lemma 1 A characteristic vector 5 k (r,p), entailed by an arc {r,p) k G 
U k \U^, where k G K is fixed, is a solution of the homogeneous linear system 



E 



x %- E 4 = °> ieI " 



ie/+(i/ fc ) jei-(u") 

Proof. Let a support Ut = {U^, k G K} be defined. For a fixed k G K 
we consider the set U^ which is, according to Theorem 2, a spanning tree 
for the network S k , and let L k p be the unique cycle of the network S k = 
(I k , U^ (J(t, P) k )i which appears after appending the arc (r, p) k G U k \U^ to 
the set f/|. 

Consider the vector x k = (x k ,, (i,j) k G U k ) of unknowns in system (81). 

Let's let a;^- = 0, (i, j) k G U k \L k p . Thus, the system (p| can be reduced 
to 

E 4" E 4 = °' <e^y» ( 9 ) 

ie/+(Lj p ) ie/-(ife p ) 

where I(L k p ) denotes all nodes in cycle L* . 

Letting x k p = 1, from the reduced system (9), we can easily define the 
values of the remaining unknowns x k j, (i,j) k G L k p \ (r, p) k : 



x% = stgn{t,]) L ^, {i,j) k G L k \ {r,p) h 



Algorithmically, after letting x k = 1, we pass from node r to node p along 
the cycle L k , consecutively setting the unknowns x k j, (i, j) k G LjL \ (r,p) k , 
to the values of signs of the corresponding arcs within the cycle L k . 

Note, the constructed solution vector x k satisfies all the rules of Definition 
4 of a characteristic vector, entailed by an arc (r, p) k , and hence S k (r, p) = x k 
is a solution of the homogeneous linear system (pi). □ 

Lemma 2. The set {S k (r, p), (r, p) k G U k \ U^} of characteristic vectors, 
where k G K is fixed, forms the basis of a solution space for the homogeneous 
system ^. 

Proof. According to Lemma 1, each characteristic vector satisfies the 
homogeneous system ^. 

By Theorem 2, for a fixed k G K, the set U^ is a spanning tree for 
the network S k = (I k ,U k ), hence \U^\ = \I k \ — 1. Thus, the number of 
characteristic vectors in the set {5 k (r, p), (r, p) k G U k \Uif\ equals \U k \Uj,\ = 
\U k \-\I k \ + l. 

Now it suffices to show that all the vectors in the set are linearly inde- 
pendent. 

Each characteristic vector S k (r, p), entailed by some arc (r, p) k G U k \U^, 
always has one and only one component, corresponding to the set U k \U^, 
that is equal to 1. It corresponds to the arc (r, p) k G U k \U^ that has entailed 
this vector. All the other components, which correspond to arcs U k \ L k , 
are equal to 0. This fact implies that any two characteristic vectors, entailed 
by different arcs, are linearly independent. □ 

Theorem 3. The general solution of system (p)]), for a fixed k G K, can 
be represented using the following form: 



4= J2 x k p sign(i,j) L "p+ lx k j- J2 x k p sign(i,j) L rp J , 

(T, P ) k GU k \U^ \ (T, P ) k £U k \U k / 

(10) 
(t,j) k eU k } x k p eR, (r, P ) k eU k \U k } 

where x k = (x k j, (i,j) k G U k ) is any partial solution of the (non-homogeneous) 
system (pi); x k p are independent variables corresponding to arcs (r, p) k G 
U k \U$. 

Proof. Let x k = (x k j, (i,j) k G U k ) be a general solution, and 
~k _ (^x k ^(i,j) k g U k ) - a partial solution, of the system (J5|. Since, by 
Lemma 2, the set {S k (r, p), (r, p) k G U k \ U^} of characteristic vectors forms 



the basis of a solution space for the homogeneous system (|8j), we can write 
the expression for x k in the following vector form: 



:r fc 



E <4 p 8\t, P )+x\ (11) 

(r,p) k £U k \U k 



as a sum of a general solution of the homogeneous system (|8j) and a partial 
solution of the non- homogeneous system (p5|); q^ G R are coefficients of the 



linear combination of characteristic vectors in (11 ). 



Rewriting (11) in the component form we obtain: 






J2 <4(r,p)+4, (i,j) k eU T ; (12) 



(T, P ) k eu k \u, k 



T 



x k =a k +x k (r,p) k eU k \U k . (13) 



From equations (13) we find a k p = x k p — x k p , ij,p) k G U k \U T and 



substitute into ( |12[ ). Finally, rewriting components of characteristic vectors 
according to (Q, we obtain the expression (10) for the general solution of the 
system ^. D 

Remark 2 In practice, for construction of a partial solution 
x k = (x k j,(i,j) k G U k ) of the system (pi), we a priori assume x k p = 0, 
(r, p) k G U k \U^ and solve the system 

/ j X ij ~ / j X ji = a ii 1^1- 

j^i+(u k ) jeir -(a*) 

Thus, formula (JTol) gets to the form: 



(T,p) fc e[/ fc \(7£ 



.<■:„ G R, (r, p) fc G £/ fc \£/£. 






Further, we will use the formula (14). 



3 Decomposition of the system 

Let Ut = {U T , k e -^} De a support of the network S for the system (ll 



We define a set C/ c == {U c C U k \U T , k G X}, |C/c| = g + |C/ | of cyclic arcs 
by selecting g + |C7o| arbitrary arcs from the sets U k \U T , k G iT. We denote 
t/jv = {^Jv, k € -K"}, t/jv = U k \(U T \J U c ), k e K - the set of remaining arcs, 
which were not included neither to the support Ut, nor to the set of cyclic 
arcs Uc- 

8 



Let's substitute the general solution (14) of the system (|5|, for each 
ke K, into Q: 

Z_v Z_v A *i Xi J = Z_v Z_v y x *i = 



(i,j)ec/ keK(i,j) 



k&K (i,j) k £U k 



E E %4 + E E *S< 

fcetf {i,j) k &U k 



E E ^ 

k£K (ij)k £U k 



k&K ( T ,p) k £U k \U k 



E x k Tp sign(i,j) L -p +x k 



ij 



(r,p) k &U k \U k 



+ 



+ E E X r P p X r P = «P' P =1 '^ 



(15) 



We change the summing order in (15): 



E E < E ^n(u) L - + E E A "/4-+ 



k£K {r,p) k £U k \U k (i,j) k &U k 



k&K (ij)k &U k 



+ E E A S< = «P' P =1 >9- 

k&K ( T ,p) k eU k \U k 



(16) 



In equations ( 16 ) we group the variables, corresponding to the sets U k \U^, 
keK: 



E E 

k&K ( T ,p) k eU k \U k 



X 



Tp 



A 5 + E A 5^n(*, 

(i,j) k eu k 

kp~k 



J) L - 



keK (i,j) k eu k 



Definition 5 We call the number 

R p( L T P 



E X ij si 9 n (hJ) 

(i,J) k £L k 



(17) 



;is) 



the determinant of the cycle L Tp , entailed by an arc (r, p) G U \U^, with 
respect to the equation with the number p of the system 



Let's denote 



^ = «--E E \Y*«> /' = i-</. 



P / 4 / J ' 'V ~IJ 

k£K (ij) k eu k 



(19) 



9 



The equations (fl7l), according to formulae (18), (19), get to the form 



fcG^(r,p) fe e(7 fe \!7*; 



(20) 



In (120]) we group the variables, corresponding to the sets Uq, k € K: 



E E R P (L k Tp )x k Tp = A?-J2 E Rp(L k Tp )x k Tp ,p = l,q. (21) 



Now, we apply the similar considerations to the system d3l). Note, that 



3) can be regarded as a particular case of the system (2) with \? equal to 
Oorl. 

Let us substitute the general solution (14) of the system ([5]), for each 
k e K , into @: 



E 4= E 4+ E 



X 



U 



k£K (i,j) keK (i,j), keK (i,j) 

Jk 



(i,j) k £U k (i,j) k £U k \U k 



£ 



(i,j) k £U k 



E x k rp sign(i,j) L -p + x k 



i.i 



(T,p) k eU k \U k 



+ 



+ E 4 = %> {hj)£u . 



k&Ko(i,j), 

(i,j) k eu k \u k 



(22) 



Now, after changing the summing order and grouping the variables, cor- 
responding to the sets U k \U^, k e K (i,j), (i,j) G Uq in (22), we obtain 



E 



xi 



x k Tp sign(i,j) L ^ = Zij- E "v 

kGKo(i,j), k£Ko(i,j), 

{T, P ) k £U k \u k (i,j) k eu k 



{i,j)eU . (23) 



On this step let us introduce the following notation: 



Sij(L k Tp ) 



sign(i,j) L -p,k e K (i,j) 
0,k<?K {i,j) 



(t,j)eU ,(r, P ) k eU k \U k ,keK. 

(24) 
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Thus, equations p3J) get to the form 



E E s n( L r P )4 P = ^ij, ihj)eu , 

kEK ( T ,p) k eU k \U k 



where 



.4 



ij 



z 'j 



E 

(i,j) k &U k 



x k - 



ihj) eu . 



(25) 



(26) 



In (25) we group the variables, corresponding to the sets U c , k G K: 

E E wo, 

k&K ( TiP )fe e c/fe 



)x k 



^-E E ^< L 






ihj) e^o- 



(27) 



Finally, let us rewrite equations (21 ) and (27) in the matrix form. For this 
purpose, we introduce arbitrary numberings of arcs within the sets Uq and Uc- 
Thus, £ = £(i,j) is a number of an arc (i,j) G £/o, £ G {1, 2, . . . , \Uq\}; and t = 
t(r, p) k is a number of a cyclic arc (r, p) k G £/£, A; G K, t G {1,2, ... , |£/c|}- 
In other words, we number the equations of the system ([3]), or (27), and 
the variables, corresponding to the set Uc- Note, the numbering of cyclic 
arcs is equivalent to the numbering of the set {L k p , (r, p) k G U c , k G K} of 
cycles, entailed by arcs (r, p) k G U c , with respect to spanning trees [/* of 
the networks S k . 



Now equations (21) and (27) can be regarded as following: 

Dx c = (3, 

l^q,t(r,p) k 



(2* 



where D 



Dl \D 1 



D ,. ,,, - (R p (L k p ),p 



l,\Uc\) - sub matrix 



of the size q x \U C \, D 2 = (<^(L* ),£(i, j) = 1, \U \,t(r,p) k = 1,\U C \) 



submatrix of the size E/q x \Uc\, xc 



(x 



t p' j 



(r,p) k eU c ,keK)- vector of 



unknowns with components ordered according to the numbering t = t(r, p) 



The right-hand side of (28) has the form: 






(29) 



where j3 p = A p 



E E R P ( L 

keK ( TjP )fe e (7fe 



Tp) X TplP 



1>?> 






TpJ^Tpl 



ihj) e ^o- 



n 



From ( 28 ) , in case of non-singularity of the matrix D, we find the unknown 



variables xc, corresponding to the set Uc of cyclic arcs: 

x c = D~ 1 p. (30) 

Remark 3. Generally, because of an arbitrary selection of arcs for the 
set Uc = {U c ,k G K}, non-singularity of the matrix D is not guaranteed. 
In the case when det D = one should re-select arcs into the set Uc and 



re-compute D,(3 for the system (28) 



Let D 1 = (vi }S ; I, s = 1, |£/c|)- We rewrite (30) in the component form: 



4 P = Yl Ut >pPp + Yl "^w^Mi * = *( r > P) fc ' ( r ' P) k £U c ,ke K. 

Thus, we have determined all the unknown variables x k = (xf-, (i,j) k G 
U k , k G K) of the system Q - (g): 

q 

X rp = / J Vt,p(3p + 2 ; U t,<l+Z(iJ)Pq+Z(i,j)) 

t = t(r,p) fc , (r,p) fc G^, fceif, (31) 



4 = ^ x* p «pn(«, j) L - + ^ + 5J-, («, JT eDj,*6 #, (32) 



x^GR,(r,p) fe G^, 
where ^ = /^ x k p sign{i 1 j) L ' r p. 

(r,p) fe e(7§ 

Note, the components of the vector x k = (x k p (i,j) k G U k ) of a partial so- 
lution of the system ^ are constructed according to the rules in the Remark 
2. 

Before we start with a simple example, let us briefly discuss the most 
important, in our opinion, aspects of the method. Although the strict esti- 
mate of complexity was left beyond the scope of the paper, one can notice 
that the described approach, if implemented on proper data structures, leads 
to efficient algorithm: the reasonable part of computations is done on small 
subsets of arcs, e.g. on 'isolated' cycles - (J7|), ( [18] ), or spanning trees - (19) 



(26). The use of the embedded network structure allows performing decom- 



position of the system and, finally, inverting the matrix D (28) of a size 
much smaller than that of the initial system pj)-(j3]). Moreover, the fact that 
the same results were obtained for each type of flow k G K, e.g. Theorem 2, 



12 



Lemmas 1 and 2, formulae (14), makes the method ready for implementation 
in parallel environment. 

However, the power of the approach is appreciated in the context of large 
problems of non-homogeneous network flow programming with (fl])-(p]) being 
the system of main constraints, where the presented ideas provide the uni- 
form technique for computing essential quantities: increment of an objective 
function, feasible directions, pseudo-flow, etc. 

Currently the authors work on the application of the obtained results for 
derivation of an optimality criterion for a broad class of non-homogeneous 
network flow programming problems. 

4 Example 

Let us consider the example (la) - (3a) of the problem ([I]) - ^ for the net- 
work S=(I,U),I= {1, 2, 3, 4, 5}, U = {(1,2), (1, 3), (2, 3), (2, 4), (3, 4), (4, 5), (5, 3)}. 
Let K = {1, 2, 3} be the set of types of flow, U 1 = {(1, 2), (1, 3), (2, 3)}, U 2 = 
U 3 = {(2,3), (2,4), (3,4), (4,5), (5,3)} - the sets of arcs carrying the flow of 
type k,k G K. We construct the networks S k = (I k , U k ), k G K (Figure [l]). 









X 12 + X 13 = ^ 














X 23 — X 12 = " 














— x 13 — x 23 = —10 












a '23 ~r x 24 — u 
2 2 2 
x 34 x 23 x 53 ~~ 

2 2 2 

x 45 x 24 x 34 — 

iy^ iy^ "I 

x 53 x 45 ~~ L 


-5 

1 












a '23 ~r x 24 — u 
3 3 3 
x 34 — x 23 — ^53 ~~ 
3 3 "i 
x 45 x 24 x 34 — 


-7 
1 












x 53 x 45 ~~ v 








It 1 


+ 3^ 


\ Xoo ~~ \ ^.Jucyn \~ £djno \~ ddjyA 


^U/ 24 


+ 2x 34 


+ x 34 






x 45 


+ 7xl 5 + xj 3 + 2xl 3 -- 


= 69 






X 12 


+ 2x 13 


+ 2x 23 
2r 2 


1 t^o^ I Qt"^ nt* 

-|- ox 23 t ox 23 x 24 

T OX 45 -|- ZX 53 X 53 


™0 

x 24 

= 58 


+ x 34 + 


x 34 



;ia) 



(2a) 



3^24 * *^24 — (^O&J 



13 



x xX 


> 


A v '•' 5 


1 
J 


^X 


---"V^-' // 




>-z--" \\ 




Jtr \ X // 








■> - 1. 





Figure 1: Union of networks S fc = (I k , U k ), k G K = {1, 2, 3} 

We choose a support of the network S 1 = (J, U) for the system (la). By 
Theorem 2 (Network Support Criterion), we build spanning trees 
U$,keK = {1,2,3}: U l T = {(U) 1 ,^) 1 }, C/| = {(2, 3) 2 , (2,4) 2 , (4, 5) 2 }, 
t/| = {(2,4) 3 ,(3,4) 3 ,(4,5) 3 }. 

Now, we compute the set {<5 fc (r, p), (r, p) fe E U k \ Utf] of characteristic 
vectors with respect to the constructed spanning tree U^, k G K = {1, 2, 3}. 



Table 1 The set of characteristic vectors with respect to the spanning tree 



XJ\ 



(ijy 


(U) 1 


(M) 1 


(2,3) 1 


*S(t,aO = <%(2,3) 


1 


-1 


1 



ul 



Table 2 The set of characteristic vectors with respect to the spanning tree 
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(U) 2 


(2,3)2 


(2,4)2 


(4,5)2 


(3,4)2 


(5,3)2 


*J(r,p) = ^(3,4) 


1 


-1 





1 





^•(r,p) = 5|(5,3) 


-1 


1 


1 





1 



Table 3 The set of characteristic vectors with respect to the spanning tree 



L/£ 



(*,j) 3 


(2,4) 3 


(3,4) 3 


(4,5f 


(2,3) 3 


(5,3) 3 


^.(r,p) = ^.(2,3) 


-1 


1 





1 





*6(r,p) = <5J(5,3) 





1 


1 





1 



Let's compute the partial solution of the system (la) for each k G K 



{1,2,3} according to the Remark 2: x 1 



lX 1 2,X 1 3,X 2 3 y l 



-6,10,0) J 



^y»^ ( ,-)-* t-f>£ rv*£ rf^ rp^ \ J- (tX f\ "1 Ci C\\ -L /-y»0 I rfJ rfJ ry& rfJ rp<J \ ■* 

x ~~ l x 23' x 24' x 45' x 34' -^53/ ~~ W) u > ■*-> u > U J 5 x ~~ l x 24' x 34' x 45' x 23' x 53/ 

(5, -7, -1,0, Of. 

3 

We form the set U C =[JU^ = {(2, 3) 1 , (3,4) 2 , (2, 3) 3 } of cyclic arcs. The 

fe=i 

3 

remaining arcs will be included into the set Un = M Un — {(5, 3) 2 , (5, 3) 3 }. 

fc=i 
Structures, representing the union of the sets U^ [J Uq, k E K — {1,2, 3} are 
shown on Figure 2. 




® 



IT 



dT 



&^ 



Figure 2: Sets U% U E/£ for networks S k ,keK = {1, 2, 3} 



Using formula (18) we compute the determinants of the cycles L k , en- 



tailed by the arcs (r^p) k E U k \U^, for each k E K = {1,2,3}, with respect 
to the equation (2a) with the number p — 1, 2 (Table 4). 

Table 4 Determinants of the cycles L^ , entailed by the arcs {r,p) k G 
f/ fc \C/|, A; e /r = {1, 2, 3} 



(r,p) k 


(2,3)! 


(3,4)2 


(5,3)2 


(2,3) a 


(5,3) 3 


Rx{L%) 





3 


-1 


7 


10 


R2(L k p ) 


1 


7 


-6 


5 


3 



15 



Now, let's compute the values 5jj(L!j. p ), (i,j) G U , (r, p) k e U k \ Uj,, k e 
K = {1,2,3} according to the formula (24) for the example (la)-(3a), U = 
{(2,4)}, 7^(2,4) = {2,3} (Table 5). 



(r, P ) k 


(2,3)! 


-vv. rpn 

(3,4) 2 


(5,3) 2 


(2,3) 3 


(5,3) 3 


M^p) 





-1 


1 


-1 






Before assembling the matrix D of the system (28), let's number the arcs 
of the set : U c = {(2, 3) 1 , (3, 4) 2 , (2, 3) 3 } : t(2, 3) 1 = 1, t(3, 4) 2 = 2, t(2, 3) 3 = 
3. The numbering within the set Uq = {(2,4)} is trivial: £(2,4) = 1. 

First, we construct the matrix D\ = (R p (L k p ), p = 1,2, t(r, p) k = 1,3) 
of the determinants of the cycles L k , entailed by the arcs (r, p) k e Uc, by 



selecting the corresponding columns from the Table 4: 



3 7 

1 7 5 



Similarly, by selecting the corresponding columns from the Table 5, we 
form the matrix D 2 = {6u(L* p ), £(2,4) = 1, t(r, p) k = i\3): 

D 2 = (0 - 1 - 1). 

Thus, joining D\ and D 2 together, we obtain the matrix of the system 



(28): 



D 




5 ) ,det£>^0. 



Let us compute the vector (3 in the right hand side of (28) using formulae 



(29): 



Pi = 
The values i? p (L 2 3 ), 



2 - 
v 53, 



X 



53 



R X {L 



x 



-A l -R x {Ll 

- A 2 - R 2 (LJ 3 )xl 3 - R 2 (Ll 
A 2i - ^24(^53)^53 _ ^24(^5 

i? p (L 3 3 ), p = 1, 2 of the determinants of the cycles 



53: 

7 53) X 53' 



J 53J X 53- 



L k p , entailed by the arcs ij,p) k e Un, as well as the values S^L^) and 
£24 (£53), are already computed and stored within the Table 4 and Table 5. 



The numbers A X ^A 2 ,A 2 ^ are evaluated using the formulae (19) and (26): 



A 1 = at-XUxlr 



A 13 x 13 



21~2 



\ zl r z -\ zi r z -\ zi r 



i21~2 



v 23" t '23 



24-^24 



M5^45 



.\31~3 . 

/v 24 x 24 



_\ 31 r 3 _A 31 r 3 — fifi 



v 34" t '34 



M5" t '45 



16 



A 1 



>12~1 . 
v 13 x 13 



\ 22 f 2 - 
x 23 x 23 



^22-2 



rv — X 1 ^ 1 — \ iz ^ 1 — \ zz ^ z —\ zz ^ z —\ zz ^ z —\ sz ^- i — \ iz ^ _\- iz Z. 
«2 A 12 x 12 A 13 x 13 A 23 x 23 A 24 x 24 A 45 x 45 A 24 x 24 A 34 x 34 A 45 x 



22~2 
45 x 45 



1 32-3 



v32-3 . 
x 34 x 34 



v32~3 



M5" t '45 



36, 



.4 



21 



^24 



r,2 



X 



7.3 



21 



X 



21 



-4. 



66 + x 2 



Thus, we have defined the vector (3 



-,:; 



10x 3 3 



36 + 6xL - 3x 



53 



-4 



x 



53 



Since the matrix D turned out to be non-singular, we can use formula 
(30) for finding the solution xc = (x£ p , (r, p) k e Uq, k G K) of the system 
28b: 

/ 1 1 l l\ 



x 



23 

2 

31 



23 



f 



V 



1 
o 
o 



66 + x\ % — IOX53 



'53 
,2 
-53 



36 + 6a4 - 3xjj 3 



-4 



x 



53 



4/ 



Finally, using formulae (31) - (32), we can define the solution of the 



system (la)-(3a) with x\ z ,x\^ being independent variables: 

27 1 



X90 — OO ZXco OJjcq. X 



23 



^53 



^53' ^34 



~2~ ~^~ 2 Xs3 ~^~ o X53,X23 



O^ rr><$ 

2 2 53 2 53 ' 



X 12 



2J 2x 53 ox 53 , 



x 1 



13 



x 



23 



x 



21 



-25 + 2x^3 


+ 8x 3 3 , 


9 
~2 


1 2 

+ 2 X 53 


5 3 
+ 2^53, 


19 

y 


1 2 
2 ^53 


5 3 
2 53 ' 



X 



45 



X 



53 



+ 1, 



X 24 



17 
-y + ^x^ -i ~> 

13 



53 



^53' 



X 



31 



Jb C 



53 



-X 



53- 



X 



45 



X 53 1' 



x 53' x 53 *= **" 
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